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Abstract 


In  this  report,  I  present  the  designs  of  two  diffractive  optical  array  generators 
used  in  an  optical  analog-to-digital  converter.  The  designs  follow  a  proce¬ 
dure  that  I  developed  for  the  design  of  general  diffractive  optical  elements. 
I  consider  both  separable  and  nonseparable  approaches  to  designing  a  5  x  5 
binary-phase  fanout,  and  discuss  indirect  and  direct  optimizations  for  the 
nonseparable  approach.  I  use  a  nonseparable  approach  with  direct  optimi¬ 
zation  to  design  an  8-phase-level  7x7  array  generator. 
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1.  Introduction 


Diffractive  optical  elements  (DOEs)  are  thin  optical  elements  with  func¬ 
tions  that  derive  from  diffraction,  as  opposed  to  reflection  or  refraction. 
Diffractive  optics  both  supplement  and  complement  refractive  optics.  That 
is,  although  a  DOE  can  be  used  to  replace  refractive  components,  it  can 
also  be  used  in  combination  with  refractive  optics  to  enhance  the  overall 
performance  of  an  optical  system.  Most  important  is  that  DOEs  can  be 
designed  and  fabricated  to  implement  functions  that  refractive  and  reflec¬ 
tive  components  cannot  implement. 

The  widespread  application  of  DOEs  requires  an  understanding  of 
diffractive  design.  Although  archived  publications  on  this  topic  are  avail¬ 
able  [1,2],  they  typically  concentrate  on  theory  and  provide  few  examples. 
The  objective  of  this  report  is  to  highlight,  through  example,  the  important 
issues  in  diffractive  design.  The  two  DOEs  I  discuss  here  were  produced 
for  an  optical  analog-to-digital  (A/D)  converter  designed  by  MAJ  Barry 
Shoop  of  the  U.S.  Military  Academy  [3]. 

Section  2  includes  a  brief  review  of  diffractive  design,  and  sections  3  and  4 
show  design  examples. 
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2,  Diffractive  Design 

The  design  of  a  DOE  incorporates  three  primary  steps:  (p  understanding 
the  physics  of  the  design  problem  (analysis),  (2)  translating  that  physical 
understanding  into  mathematics  and  defining  an  appropriate  optimization 
problem  (sjmthesis),  and  (3)  executing  the  design  and  fabricating  the  ele¬ 
ment  (implementation).  Following  is  a  further  breakdown  of  each  step. 

2.1  Analysis 

Step  1:  Understand  the  physics  of  image  formation  using  the  DOE. 

Step  2:  Understand  the  fabrication  of  the  DOE  from  data  generated  by 
computer. 

By  ^^understand'^  I  mean  to  imply  the  selection  of  a  model.  The  coherent 
Fourier  transformer  shown  in  figure  1  is  a  valid  representation  of  the 
optics  for  the  two  examples  that  I  consider  here.  Under  the  assumptions  of 
scalar  diffraction  theory,  the  DOE  is  a  physically  thin  element  that  modu¬ 
lates  the  complex  wave  amplitude  of  an  incident  wavefield  according  to 

P(u,  v)  =  A{u,  v)  exp  [jO  (u,  v)]  .  (1) 

If  the  sotuce  illumination  is  a  spatially  coherent,  quasi-monochromatic 
plane  wave  with  wavelength  A,  a  Fourier  transform  relatioriship  exists 
between  the  diffractive  element  P(u,  v)  in  the  pupil  plane  and  its  response 
p{x,  y)  in  the  image  plane  [4], 

r  r  Stt 

p{x,y)  =  J  J  P{u,v)exp  jjj{ux  +  vy)  dudv  *  (2) 

where /is  the  focal  length  of  the  lens. 

Assuming  that  a  multistep  etch  procedure  is  used  to  fabricate  the  DOEs 
[5],  if  the  fabrication  process  from  computer-generated  data  to  diffractive 
element  is  linear  and  introduces  no  errors,  I  can  model  the  DOE  as  a 
multilevel  quantized  phase  element  (i.e.,  A{u^v)  =  1  and  6{u^v)  - 
27r£/L,  ^  =  [0;L-1]). 

Figure  1.  Coherent 
optical  Fourier 
transformer. 


have  ignored  complex  scale  factors. 
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■pie  pOEs  that  I  considered  each  produce  an  array  of  point  sources  when 
illuminated  by  a  single  point  source  and  are,  therefore,  referred  to  as  array 
generators.  I  represent  the  desired  array  by  q{x,  y)., 

N 

Q(^,  y)  =  yn)  >  (3) 

n=l 

where  N  is  the  total  number  of  spots  in  the  array,  (xn,  y„)  are  their  spot 
locations,  and  the  phasor  qn  is  their  spot  magnitude  and  phase.  If  the  spot 
locations  (x„,  y„)  are  rational,  the  Fourier  transform  Q{u,  v)  of  q{x,  y)  is 
periodic  and,  therefore,  P(u,  u)  must  also  be  periodic, 

P{u,v)  =  P(u,v)  *  *  repl(M,'y),  (4a) 

where  p{u,  w)  is  a  single  period  representation  of  the  array  generator,  with 
unity  extent  along  both  axes.  The  two-dimensional  convolution  (repre¬ 
sented  by  *  *)  of  P[u,  v)  with  the  replicating  function  repl(u,  v)  tiles  the 
Fourier  plane  (see  fig.  2).  As  a  consequence,  the  points  in  image  space  at 
which  the  generated  array  has  value  are  determined  by  the  lattice  array 
samp(a:,y). 


P{x,  y)  =  p{x,  y)  samp(a;,  y).  (4b) 

The  functions  repl(«,  v)  and  samp(rc,  y)  form  a  Fourier  transform  pair. 


2.2  Synthesis 

Step  3:  Define  the  design  metric. 

Step  4:  Define  the  optimization  problem. 

The  fidelity  of  an  array  generator  is  based  on  its  ability  to  generate  the 
desired  spot  array  with  high  diffraction  efficiency  and  minimum  intensity 
error.  Diffraction  efficiency  is  the  ratio  of  the  energy  diffracted  into  the 
desired  signal  to  the  total  energy  in  the  image  plane. 


/  /  \p{x,y)\^dxdy  f  f  \p(x,y)\^dxdy 

r,  =  _  _  (x,y)ex 

/  oo  oo  5o  oo  ” 

I  I  \p{x,y)\^dxdy  J  f  \P{u,v)\'^dudv 

-OO  ~oo  _oo  ^OO 

/  J  \p{x-,y)?dxdy  =Y^\Pn?- 

{x,y)ex 

Intensity  error  is 

eint=  j  j  \\p{x,y)f  -  a\q{x,y)'^f  dxdy  = 


(6) 
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My  use  of  the  scale  factor  a  indicates  that  the  absolute  intensity  of  the 
response  is  not  critical;  rather,  the  response  p[x,y)  should  exhibit  the 
shape  described  by  q{x,  y).  This  is  defined  as  an  application  of  scale  free¬ 
dom  [6].  However,  a  scale  factor  Omin  exists  that  minimizes  eint: 


f  I  \p{x,y)\  \Q{x^y)\  dxdy  ^  \pn?\qn? 

(x,y)eX _ _  _  n=l _ (7) 


f  f  \Q{^,y)\  dxdy 

(x,t/)€X 


N 

n=l 


The  scale  a^in  yields  the  minimum  fluctuation  from  the  desired  intensity 
by  the  generated  intensity,  which  is  represented  in  figure  3. 

The  physical  extent  of  the  detector  in  the  image  plane  and  the  fact  that  it 
responds  to  intensity  influence  the  definition  of  the  error.  I  am  concerned 
only  with  the  intensity  of  the  response  within  the  signal  window  x;  the 
phase  within  the  signal  window  and  the  total  response  outside  of  the  sig 
nal  window  are  unimportant.  Thus,  I  am  also  using  phase-freedom  and 
complex-wave  amplitude  in  my  designs  [6].  However,  a  tradeoff  exists 


Figure  2.  Transform 
relationship  between 
Fourier  replicating 
function  and  lattice 
of  image  sampling 
points  that  it 
generates:  (a)  tiling  of 
Fourier  plane  with 
P  (u,  v)  by  repl  (u, v); 
(b)  lattice  image 
sample  points 
samp  (a:,  y)  by 
repl  (wjv). 


(a) 


(b) 


X 


Figure  3.  Significance 
of  scale  factor  amin  in 
determining  intensity 
error. 
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between  diffraction  efficiency  and  error  (i.e.,  minimum  error  can  be 
achieved  by  reducing  diffiraction  efficiency).  Likewise,  high  diffraction  effi¬ 
ciency  can  be  achieved  by  allowing  for  large  levels  of  error.  One  deter¬ 
mines  these  measures  after  the  design  is  complete. 

For  optimal  design,  it  is  necessary  to  define  a  metric  that  accoimts  for  this 
tradeoff  before  beginning  the  work,  and  to  use  it  within  an  optimization 
routine.  I  have  used 


^ub 


N 

n=l 


\Pn\  ^ub  i^nl 


(8) 


where  =  r].ub/N  and  rjub  is  the  diffraction  efficiency  upper  bound  [7], 


?7jx6  — 


1/2  1/2 

/  /  \Q{u,v)\dudv 

-1/2  -1/2 


n  2 


1/2  1/2 

/  /  \Q{u,v)'^ dudv 

-1/2  -1/2 


(9) 


Thus,  the  intensity  distribution  OLub\(l{,^ty)\^  represents  the  maximum 
array  intensity  that  a  phase-only  E)OE  can  generate.  Although  t]  =  100  per¬ 
cent  also  represents  an  upper  bound  on  diffraction  efficiency,  it  is  too  gen¬ 
eral  to  be  effective;  provides  more  specific  information  about  the  par¬ 
ticular  design  of  interest.  Because  Oiub\q{x,  y)\‘^  represents  a  fixed,  upper 
limit  in  the  space  of  solutions,  the  metric  is  an  absolute  distance 
measurement. 


2.3  Implementation 

Step  5:  Select  a  design  algorithm  and  performing  the  optimization. 

Step  6:  Fabricate  the  DOE  from  the  computer-generated  data. 

Step  7:  Test  the  DOE. 

Steps  5  to  7  are  straightforward — select  the  design,  fabricate  tiie  DOE,  and 
test  its  performance  in  the  system — ^but  require  many  practical  decisions 
and,  most  importantly;  a  considerable  expenditure  of  resources. 

One  of  these  decisions  is  determining  the  physical  dimensions  and  mini¬ 
mum  features  of  the  DOE  before  fabrication.  The  transform  relationship 
between  P(«,  v)  and  p{x,  y)  assumes  that  P(m,  v)  has  unit  extent  along 
both  axes.  The  spots  of  p{x,  yf  therefore,  are  located  at  integer  multiples  of 
A/-  For  an  arbitrary  spot  spacing  d  along  one  axis,  the  unit  period  width 
W  is 


w  =  Xf/d. 


(10) 
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If  the  total  length  of  the  DOE  is  D  in  one  dimension,  the  number  of  replicas 
M  of  the  unit  period  is 

M  =  DJW.  (11) 

Further,  if  the  minimum  feature  is  A,  the  number  of  addressable  features 
L  within  the  unit  period  is 

L  =  W/A.  (12) 

The  total  number  of  addressable  features  in  the  element  is  its  space- 
bandwidth  product  (SBWP)  S, 

S  ^ML  =  D/A.  (13) 

For  example,  for  a  1-cm,  one-dimensional  element  written  with  1-pm  fea¬ 
tures,  S'  =  10^.  Due  to  the  fact  that,  at  a  minimum,  it  takes  two  features  to 
code  a  single  value  in  the  desired  response,  the  SBWP  has  a  direct  imjpact 
on  design-  Therefore,  the  maximum  number  of  points  that  can  be  specified 
in  q{x,y)  is  S/2. 
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Design  Examples 

3.1  Example  1:  (0, 7r)-Binary-Phase  5x5  Fanout 

The  first  design  example  is  a  (0, 7r)-binary-phase  5x5  fanout  (a  ^^fanout^^ 
implies  that  all  spots  have  equal  intensity).  The  illumination  wavelength 
A  =  850  nm  and  the  focal  length  of  the  lens  /  =  50  mm.  The  spacing 
between  spots  along  both  axes  is  d  =  250  |xm.  The  grating  period  is,  there¬ 
fore,  170  X 170  pm2.  The  element  size  D^.  x  Dj,  =  5  x  5  mm^  which  impUes 
that  the  unit  cell  is  repeated  30  x  30  times. 

For  this  example,  two  approaches  to  design  are  possible;  separable  and 
nonseparable.  A  two-dimensional  function  is  separable  if  its  behavior 
along  one  axis  is  independent  of  its  position  on  the  other  axis.  In  a  sepa¬ 
rable  design,  one  considers  the  desired  performance  along  the  horizontal 
and  vertical  axes  separately,  and  crosses  the  one-dimensional  results  to 
construct  a  two-dimensional  array  generator.  A  nonseparable  design  con¬ 
siders  the  two-dimensional  performance  as  a  whole. 

The  design  requires  determination  of  the  upper  bound  on  the  diffraction 
efficiency.  Except  for  the  constraints  imposed  on  array  phase  by  a  binary- 
phase  array  generator,  which  I  discuss  below,  determination  of  the  diffrac¬ 
tion  efficiency  upper  bound  is  independent  of  any  fabrication  consider¬ 
ations.  The  upper  bound  represents  the  theoretical  upper  limit  on  the 
diffraction  efficiency  with  which  the  desired  array  can  be  generated. 

The  expression  for  the  diffraction  efficiency  upper  bound  requires  com¬ 
plete  specification  of  the  desired  array  q{x,  y)  in  terms  of  magnitude  and 
phase.  The  design,  however,  considers  only  the  intensity  |?(x,  y)\^,  and  so 
only  the  magnitude  y)|  is  specified.  This  limitation  is,  however,  to 
one  s  advantage  in  terms  of  high  diffraction  efficiency  design,  since  the 
freedom  to  specify  phase  allows  a  choice  of  the  phase  that  generates  the 
highest  diffraction  efficiency. 

To  determine  this  phase,  the  expression  for  the  diffraction  efficiency  upper 
bound  with  respect  to  the  array  phase  [7,8]  must  be  maximized.  However, 
a  binary-phase  DOE  generates  a  response  that  is  Hermitian  (i.e., 
y)  P  (  X,  j/));  therefore,  the  phase  of  q{x,  y')  is  not  completely 
arbitrary,  but  must  be  odd.  The  array  phase  6ub{x,  y), 


y) 

X  =  —2 

o 

II 

1— 1 

1 

II 

x  =  1 

X  =  2 

y  =  -2 

-1.3915 

1.2433  -1.2749 

1.6429 

0.7962 

y  =  -l 

-0.1648 

1.6689  2.3327 

-0.0511 

-0.3013 

y=  0 

1.5014 

-2.5544  -3.1416 

2.5544 

-1.5014 

y=  1 

0.3013 

0.0511  -2.3327 

-1.6689 

0.1648 

y=  2 

-0.7962 

-1.6429  1.2749 

-1.2433 

1.3915 
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yields  an  upper  bound  of  78.88  percent;  thus,  oiufe  =  31.5  x  10  •.  ' 

valued  Fourier  transform  Qub{u,  v)  of  <lub{^i  v)  —  Vn  u6(.^)  y)l 

is  represented  in  figure  4. 


3.1.1  Separable  Design 

The  nature  of  the  5  x  5  fanout  allows  construction  of  a  two-dimensional 
array  generator  from  a  single  one-dimensional,  5-spot  fanout.  Further,  the 
array  generator's  binary-phase  and  small  number  of  spots  allow  use  of 
so-called  Dammann  design  techniques  [2,9]. 

To  achieve  a  uniform-intensity  spot  array,  Dammann  recognized  that  a 
one-dimensional,  (0,  ;r)-binary-phase  array  generator  can  be  represented 
completely  by  the  locations  Zk  at  which  its  phase  transitions  occur.  Con- 
sider  figure  5  and  the  representation  of  P(«)  in  terms  of  a  binary  ampU- 
tude  function  E(t4), 

P{u)  =  2E(«)  -  rect(M),  (14a) 


where 


K-l 

T,{u)  =  ^  rect 

Jb=l 
k  odd 


u  —  {zk+i  +  Zk)  /2 


Figure  4. 

Representation  of 
Qub(u,v)  for  a  5  X  5 
fanout  (two  periods 
are  shown). 


(14b) 


Figure  5.  Binary 
amplitude 
distribution  S(u) 
used  to  describe  a 
one-dimensional 
binary-phase  array 
generator. 


1  r- 


1/2  zi  22  2:3  24 


25 


^K-2 
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If  repl(u)  =  comb(M),  where  comb(w)  =  §  6{u-  n),  then,  by  virtue 

of  the  self-transforming  properties  of  the  comb-function,  the  generated 
array  is 


00 

[2o-(«)  -  sinc(n)]  6{x-n) , 

n=— oo 


(15a) 


where 

K-l 

=  E  (^fe+i  -  Zk)  sine  [{zk+i  -  Zk)  n]  exp  [jtt  (^^+1  +  Zk)  n] .  ,, 

*==1  (15b) 

fc  odd 

The  comb-function  replication  generates  a  spot  array  that  assumes  values 
at  integer  locations  in  the  image  plane  and,  from  symmetry,  generates  an 
array  that  has  an  odd  number  of  spots.  To  generate  an  even  number  of 
spots,  it  is  necessary  to  use  an  alternative  replication  [10], 

P{u)  =  P{u)  *  exp(j7r«)comb(u)  , 

which  generates 

p{x)  =  p(a;)comb(a:  -  1/2). 

The  ;r-phase  ch^ge  between  replicas  acts  as  a  spatial  carrier  to  translate  all 
spots  off  am  without  changing  the  spacing  between  sample  points.  In  con¬ 
junction  with  the  symmetry  imposed  by  the  grating,  this  implies  that  the 
array  has  an  even  number  of  spots. 

To  determine  the  phase-transition  locations,  the  generated  intensity  \p{^)\^ 
proportional  to  the  desired  intensity  1^(2:)  p,  |p(x)p  = 
r  for  X  =  [0,  iVj.  Ehie  to  the  Hermitian  symmetry  imposed  by  a 
(0,  ;r)-binary-phase  array  generator,  the  constant  a  is  real-valued.  The 
N+1  real-valued  equations  can  be  solved  for  a  and  the  phase  transitions 
Zkfk  =  [I,  K].  Because  the  system  of  equations  is  nonlinear,  there  is  no 
constraint  that  the  number  of  unknowns  must  equal  the  number  of  equa¬ 
tions— the  total  number  of  transitions  K  is  arbitrary.  However,  in  general, 
the  fewer  the  number  of  phase  transitions,  the  higher  the  diffractiorl 
efficiency. 

The  phase  transitions  for  many  one-dimensional,  Dammann  gratings 
((0,  ;c)-binary  phase  gratings)  can  be  found  in  several  references  (for 
example,  [11,12]).  The  phase  transitions  within  a  unit  cell  of  length  unity 
that  generate  a  5-spot  fanout  with  high  diffraction  efficiency  are  shown  in 
table  1.  The  value  of  o;  for  the  one-dimensional  problem  is  0.1548,  the 
square  of  which  (24.0  x  IQ-^)  is  the  value  for  the  two-dimensional  problem. 
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Table  1.  Phase  transitions  in  a  Dammann  grating  to  produce  a  5-spot  array  generator. 


Unit  cell 

End 

point 

Zl 

Z2 

^3 

End 

point 

Unit  interval 

-0.50000 

-0.36766 

-0.01930 

0.01930 

0.36766 

0.50000 

170-p,m  period  (|im) 

0.00 

22.50 

81.72 

88.28 

147.50 

170.00 

Quantized  ITO-lim 
period  (p,m) 

0.0 

22.0 

81.4 

88.0 

147.4 

170.5 

Quantized  unit 
interval 

-0.50000 

-0.37097 

-0.02258 

0.01613 

0.36452 

0.50000 

Fabrication  requires  the  conversion  of  the  normalized  locations  of  phase 
transition  to  distances  through  multiplication  by  the  grating  period  and 
addition  of  an  appropriate  shift;  these  values  are  also  indicated  in  the  table. 
Unfortunately,  the  fabrication  process  does  not  have  uriimited  resolution, 
and  distances  must  be  quantized  according  to  the  minimum  placement 
resolution  and  minimum  features  of  the  fabrication  technology.  Typical 
values  for  fabrication  technology  available  today  are  O.l-pm  placement 
resolution  and  1-pm  minimum  feature.  Table  1  includes  the  phase  transi¬ 
tions  for  a  minimum  feature  of  1.1  pm.  The  impact  on  fabrication  param¬ 
eters  is  shown  in  table  2. 

The  separable  two-dimensional  pattern  generated  from  the  one¬ 
dimensional  structure  is  represented  in  figure  6.  The  fidelity  of  the 
generator's  response  is  summarized  in  table  3.  The  generator's  intensity 
response,  as  well  as  that  of  all  the  other  designs,  is  listed  in  appendix  A. 
For  a  fanout  element,  in  which  case  l9n|  =  1/  equation  (7)  reduces  to 


^ave  ^min  ^  j  \P'^\  * 

■'''  n=l 


Further,  I  define  e^ms  as 


The  ratio  erms  /^ave  is,  in  essence,  a  per-spot  measure  of  the  noise  fluctua¬ 
tions  (i.e.,  a  noise-to-signal  ratio  per  spot).  The  tradeoff  between  diffraction 
efficiency  and  error  is  clearly  evident  in  table  3. 


3.1.2  Nonseparable  Design 

In  the  previous  section,  separable  design  of  the  array  generator  used  ana¬ 
lytic  techniques  to  determine  the  locations  of  phase  transition;  however,  in 
general,  die  desired  design  metric  must  be  optimized  explicitly.  This  is  the 
approach  that  I  take  in  this  section.  The  optimization  problem  can  be 
stated  in  the  following  way. 

Problem  1:  Minimize  Bub  over  P(u,v)  =  A(u,u)  exp[j0(«,u)]  subject 
to  A(UfV)  =  1  and  d(u,v)  =  {27r£/ =  li2, ...,L—  1}. 
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Table  2.  Summary  of 
fabrication  param¬ 
eters  for  5  X  5-spot 
array  generators.*^ 

Design 

A  =  850  nm, 

/  =  50  mm 

A. 

(lam) 

Lx 

(gm) 

JV4 

(mm) 

4 

(gm) 

Desired 

— 

— 

170.00 

— 

5.000 

250.00 

Separable 

1.10 

155 

170.50 

30 

5.115 

249.27 

Nonseparable 

1.33 

128 

170.24 

30 

5.107 

249.65 

*Due  to  square  nature  of  array  generator,  features  in  only  one  dimension  are 
indicated. 


Figure  6. 

Representation  of 
Dammann-designed 
separable  array 
generator  (two 
replications  of  the 
unit  cell  are  shown). 


IB  I 


Table  3.  Summary  of 
response  characteris¬ 
tics  for  5  X  5-spot 
array  generators. 

Array  generator 

V 

{%) 

‘'ave 

(xl0~3) 

^rms 

(xlO-3) 

^rms/^ave 

(%) 

Order 

(computing  time) 

Upper  bound 

78.88 

31.5 

— 

— 

Minutes 

Dammann  grating 

59.89 

24.0 

0.06 

0.24 

Negligible 

Quantized  design 

80.58 

32.2 

9.60 

29.65 

Negligible 

Iterated  design 

69-67 

27,9 

1.70 

5.98 

Minutes 

Annealed  design 

76.34 

30.5 

1.40 

4.46 

Hours 

Two  basic  approaches  to  solving  problem  1  have  emerged  in  the  DOE 
literature — one  that  is  direct  and  one  that  is  indirect  [2],  In  an  indirect 
design,  the  imconstrained  design  problem  must  be  solved. 

Problem  2:  Minimize  e^b  over  P(u,  v)  • 

The  solution  of  problem  2  Qubi'^',  v)  is  then  mapped  onto  the  fabrication 
constraints.  Mapping  can  be  as  simple  as  quantization,  or  it  can  be  an  opti¬ 
mal  routine.  In  a  direct  design,  the  constraints  of  fabrication  must  be 
applied  directly,  as  in  problem  1. 
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Due  to  fabrication  constraints,  the  ordinate  values  (u,  v)  are  not  continu¬ 
ous,  but  discrete  (i.e.,  (u,  v)  =  (fu  A,  i^A),  where  A  is  the  minimum  fea¬ 
ture  size  and  and  iy  are  integers.  To  accommodate  the  use  of  the  fast 
Fourier  tranisform  in  the  solution  of  problems  1  and  2, 1  divided  the  unit 
cell  into  128  minimum  features,  which  correspond  to  A  =  1.33  |tm.  The  fab¬ 
rication  parameters  for  nonseparable  design  are  summarized  in  table  2. 1 
solve  problems  1  and  2  in  the  following  sections. 

3.1.3  Indirect  Design 

Indirect  design  of  the  array  generator  requires  the  solution  of  problem  2, 
which,  as  I  stated,  is  tire  complex  fimction  Qy},{u^  v).  The  second  step  in 
the  design  is  the  application  of  the  fabrication  constraints,  which,  for  this 
problem,  implies  that  the  DOE  can  assume  only  two  values:  -1  and  +1. 
This  constraint  can  be  applied  by  simple  quantization,  which  yields  the 
array  generator  represented  m  figure  7.  Its  performance,  which  is  notably 
poor,  is  summarized  in  table  3.  It  should  not  be  concluded  from  this  sum¬ 
mary,  however,  that  indirect  methods  always  produce  poor  results.  In  this 
case,  the  poor  results  were  a  consequence  of  the  parameters  particular  to 
this  design. 

3.1.4  Direct  Design 

Two  general  nonlinear  optimization  algorithms  capable  of  solving  prob¬ 
lem  1  are  represented  in  figure  8  [13].  The  operator  T  reflects  the  imder- 
standing  (from  step  1)  of  the  influence  that  the  DOE  has  on  the  image  it 
produces,  and  the  operator  E  is  the  design  metric. 

Fienup  [13]  refers  to  the  algorithms  in  figures  8(a)  and  (b)  as  error-reduc¬ 
tion  and  input-output;  however,  to  indicate  the  nature  of  data  flow,  I  prefer 
the  nomenclature  bidirectional  and  unidirectional.  Note  that  use  of  a 
bidirectional  algorithm  implies  not  only  an  understanding  of  the  operator 
T,  but  also  an  imderstanding  of  its  inverse,  T“^  (i.e.,  an  vmderstanding  of 
how  variations  in  the  response  affect  the  DOE).  Under  my  assumptions, 
the  system  T  is  modeled  by  a  Fourier  transform,  which  is  easily  inverted. 
Therefore,  figure  8(a)  represents  the  iterative  Fourier  transform  algorithm 
(IFTA)  [6,13,14]. 

Figure  7. 

Representation  of 
indirect  design  of  a 
nonseparable  array 
generator  (two 
replications  of  the 
unit  cell  are  shown). 
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Figure  8.  Optimi¬ 
zation  algorithms 
used  for  design: 

(a)  bidirectional; 

(b)  unidirectional. 


(a) 


Although  the  IFTA  is  simple  in  its  implementation,  it  is  prone  to  stagnate 
in  local  minima.  To  overcome  this  problem,  effective  use  of  the  IFTA 
requires  the  exploitation  of  the  design  freedoms  so  that  the  algorithm  can 
modify  the  complex  transmission  values  of  the  DOE,  while  maintaining 
the  desired  performance  objective  m  the  image  plane  [6]. 

By  contrast,  the  philosophy  behind  unidirectional  algorithms  is  that,  if  the 
DOE  can  be  characterized  by  a  finite  set  of  quantized  parameters  (e.g., 
phase  levels),  a  finite  (albeit  large)  number  of  permutations  of  these 
parameters  exists.  The  objective  of  the  algorithm  is  to  find  the  permutation 
that  achieves  the  optimization  in  an  efficient  manner.  Two  examples  of  uni¬ 
directional  algorithms  are  gradient  descent  methods  and  annealing 
algorithms  [15-17].  In  these  algorithms,  the  present  DOE  is  a  function  of 
the  previous  DOE  and  the  performance  generated  by  it.  When  the  model 
of  ttie  optical  system  cannot  be  inverted,  \midirectional  algorithms  must  be 
used. 

The  results  of  the  direct  solution  of  problem  1,  using  both  the  IFTA  and 
simulated  annealing,  are  represented  in  figure  9  and  summarized  in  table 
3.  Although  there  is  a  notable  improvement  m  performance  over  quantiza¬ 
tion,  it  is  again  impossible  to  make  general  statements  about  the  effective¬ 
ness  of  simulated  annealing  over  EFTA  from  this  example.  The  diffraction 
efficiency  of  the  IFTA-designed  array  generator  is  lower  than  that  of  the 
annealed  array  generator,  but  their  performances  are  otherwise  compa¬ 
rable.  It  is  true,  however,  that  simulated  annealing  requires  more  comput¬ 
ing  time  than  the  IFTA. 
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Figure  9. 

Representation  of 
direct  design  of 
nonseparable  array 
generators:  (a)  IFTA 
design;  (b)  annealed 
design  (two  replica¬ 
tions  of  the  unit  cells 
are  shown). 


3.2  Example  2: 8-Phase-Level  7x7  Array  Generator 


The  second  design  example  is  for  an  8-phase-level  7x7  array  generator 
that  has  the  nonuniform  intensity  distribution  v)  \  indicated  below. 


X  =  —3 

X  =  —2 

a:  =  — 1 

X  =  0 

X  =  1 

x  =  2 

X  =  3 

2/  =  -3 

0.0003 

0.0019 

0.0051 

0.0068 

0.0051 

0.0019 

0.0003 

y  =  -2 

0.0019 

0.0103 

0.0248 

0.0328 

0.0248 

0.0103 

0.0019 

y  =  -l 

0.0051 

0.0248 

0.0583 

0.0766 

0.0583 

0.0248 

0.0051 

y  =  0 

0.0068 

0.0328 

0.0766 

0.0000 

0.0766 

0.0328 

0.0068 

y=l 

0.0051 

0.0248 

0.0583 

0.0766 

0.0583 

0.0248 

0.0051 

y  =  2 

0.0019 

0.0103 

0.0248 

0.0328 

0.0248 

0.0103 

0.0019 

y  =  3 

0.0003 

0.0019 

0.0051 

0.0068 

0.0051 

0.0019 

0.0003 

Unlike  the  5x5  fanout,  the  spot  spacing  along  the  horizontal  and  vertical 
axes  is  unequal:  dxX  dy  =  160  x  80  |jm^.  The  wavelength  A  =  850  run  and 
the  focal  length  j  =  50  mm;  therefore,  the  grating  period  Wx  x 
265.625  X  531.25  pm^.  The  size  of  the  DOE  was  fixed  at  1  x  1  cm^,  and  so, 
for  ease  of  fabrication  and  design,  I  selected  the  minimum  feature 
Ax  y.  Ay  =  2.1  X  4.2  pm^. 

I  designed  a  nonseparable  DOE  using  a  direct  approach  that  I  imple¬ 
mented  using  the  IFTA.  I  selected  the  number  of  features  in  a  unit  cell  to  be 
128  X  128,  which  yields  a  grating  period  of  268.8  x  537.6  pm^.  The  unit 
period  was  replicated  37  x  18  times  to  produce  an  element  9.9456  x  9.6768 
mm^.  Given  these  parameters,  the  spot  spaciag  is  158.11  x  79.06  pm^.  (Note 
that  a  separable  design  is  not  possible  because  the  DOE  has  eight  phase 
levels — ^when  crossed,  an  8-level,  one-dimensional  DOE  does  not  5deld  an 
8-level,  two-dimensional  DOE.) 


14 


A  diffraction  efficiency  upper  bound  of  91.15  percent  is  generated  by  the 
phase 


y) 

X  =  — 3 

X  —  —2  X  =  —1  X  =  0  X  =  1 

X  =  2 

X  =  3 

y  =  -3 

1.2566 

0.2038  -0.1017  -2.9513  -0.1017 

0.2038 

1.2566 

y  =  -2 

1.1104 

0.5581  -0.7317  -2.7790  -0.7317 

0.5581 

1.1104 

y  =  -1 

0.6649 

-0.2590  -1.8587  2.5672  -1.8587 

-0.2590 

0.6649 

y  =  0 

-2.9513 

-2.7790  2.5672  0.0000  2.5672 

-2.7790 

-2.9513 

y  =  1 

0.6649 

-0.2590  -1.8587  2.5672  -1.8587 

-0.2590 

0.6649 

y=2 

1.1104 

0.5581  -0.7317  -2.7790  -0.7317 

0.5581 

1.1104 

y=3 

1.2566 

0.2038  -0.1017  -2.9513  -0.1017 

0.2038 

1.2566 

The  phase  assumes  an  8-level  phase-only  DOE;  however,  no  symmetry 
constraints  were  imposed  on  the  problem — the  S5ntnmetry  exhibited  in  the 
phase  results  from  the  s5TTimetry  in  the  array  intensity.  The  magnitude  and 
phase  of  the  upper  bound  Qubl^,  v)  in  the  spatial  frequency  domain  are 
represented  in  figures  10(a)  and  (b).  An  alternate  representation  is  the  com¬ 
plex  plane  representation  in  figure  11.  The  IFTA  must  take  these  data 
points  and  map  them  to  eight  equally  spaced  discrete  values  on  the  unit 
circle. 


The  phase  of  the  element  designed  by  the  IFTA  is  represented  in  figure  12. 
The  intensity  \p{x,  y)^  generated  by  this  element  follows. 


\p{x,y)\'^ 

X  =  —3 

X  =  — 2 

X  =  —1 

X  =  0 

X  =  1 

X  =  2 

X  =  3 

y  =  -3 

0.0003 

0.0016 

0.0046 

0.0062 

0.0046 

0.0016 

0.0003 

y  =  -2 

0.0016 

0.0093 

0.0235 

0.0307 

0.0235 

0.0093 

0.0016 

y  =  -1 

0.0046 

0.0234 

0.0541 

0.0698 

0.0541 

0.0234 

0.0046 

y  =  0 

0.0062 

0.0305 

0.0699 

0.0000 

0.0699 

0.0305 

0.0062 

y  =  1 

0.0046 

0.0234 

0.0541 

0.0698 

0.0541 

0.0234 

0.0046 

y  =  2 

0.0016 

0.0093 

0.0235 

0.0307 

0.0235 

0.0093 

0.0016 

y  =  3 

0.0003 

0.0016 

0.0046 

0.0062 

0.0046 

0.0016 

0.0003 

The  diffraction  efficiency  of  the  generated  array  is  91.97  percent,  which 
is  greater  than  the  upper  bound  due  to  the  presence  of  error:  Cint  =  6.90  x 
lO"^  for  amin  =  0.92.  The  noise  fluctuations  per  spot  erms  =  0.376  x  10“^. 
Therefore,  the  DOE  satisfies  tire  design  specifications  exceptionally  well. 
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Figure  10. 
Representation  of 
Qub{u,v)  for  a  7x7 
fanout  (two  periods 
in  the  spatial  fre¬ 
quency  domain  (u^v) 
are  shown  in 
(a)  |Qw6(^9v)|,and 
<b)  arg  {Qu6(it,v)}). 


Figure  11.  Complex 
plane  representation 
of  Qub{u,v), 


Figure  12.  Phase 
representation  of 
8-level  phase-only 
array  generator  (two 
replications  of  the 
unit  cell  are  shown). 


Fabrication  of  the  8-level  element  requires  the  three  binary  masks  shown  in 
figure  13.  White  indicates  the  areas  to  be  etched;  black  indicates  the  areas 
in  which  etching  should  not  take  place.  The  horizontal  scale  of  the  masks  is 
half  that  of  the  phase  function,  due  to  the  different  horizontal  and  vertical 
spacings  between  spots.  The  mask  set  was  used  by  Honeywell  to  fabricate 
the  element,  which  is  currently  being  characterized  at  the  Army  Research 
Laboratory. 


4.  Summary 

In  this  work,  I  used  two  applications  of  DOEs  in  an  optical  A/D  converter 
to  highlight  the  procedure  for  designing  diffractive  elements  using  scalar 
diffraction  theory.  The  selection  of  models  for  the  optical  system  and  fab¬ 
rication  technology,  as  well  as  the  proper  statement  of  the  design  in  terms 
of  optimization,  is  essential  to  this  procedure.  I  designed  high-efficiency, 
low-noise  elements  using  phase,  amplitude,  and  scale  design  freedoms. 


Figure  13.  Binary 
amplitude  mask  set 
necessary  to  fabricate 
an  8-level  phase-only 
array  generator 
(single  period  repre¬ 
sentation  of  mask 
used  to  etch 

(a)  TT-phase, 

(b)  7r/2-phase, 
and  (c)  ^/4-phase). 
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Appendix  A 


Appendix  A 

The  following  tables  indicate  the  intensity  responses  of  the  array  genera¬ 
tors  designed  in  section  3. 

Upper  Bound: 


X  =  —2 

X  =  —1 

a:  =  0 

X  =1 

x  =  2 

y  =  -2 

0.0315 

0.0315 

0.0315 

0.0315 

0.0315 

y^-1 

0.0315 

0.0315 

0.0315 

0.0315 

0.0315 

y  =  0 

0.0315 

0.0315 

0.0315 

0.0315 

0.0315 

y=l 

0.0315 

0.0315 

0.0315 

0.0315 

0.0315 

y  =  2 

0.0315 

0.0315 

0.0315 

0.0315 

0.0315 

Dammann  Grating: 

\p{x,y)\^ 

X  =  —2 

X  =  — 1 

a:  =  0 

X  =  1 

X  =  2 

y  =  -2 

0.0240 

0.0239 

0.0240 

0.0239 

0.0240 

y  =  -1 

0.0239 

0.0239 

0.0239 

0.0239 

0.0239 

y  =  0 

0.0240 

0.0239 

0.0240 

0.0239 

0.0240 

y=l 

0.0239 

0.0239 

0.0239 

0.0239 

0.0239 

y  =  2 

0.0240 

0.0239 

0.0240 

0.0239 

0.0240 

Quantized  Upper  Bound: 

\pix,y)\^ 

X  =  — 2 

X  =  — 1 

X  =  0 

X  =  1 

x  =  2 

y  =  -2 

0.0225 

0.0265 

0.0239 

0.0365 

0.0294 

y  =  -1 

0.0244 

0.0347 

0.0333 

0.0532 

0.0454 

y  =  0 

0.0187 

0.0395 

0.0298 

0.0395 

0.0187 

y  =  1 

0.0454 

0.0532 

0.0333 

0.0347 

0.0244 

y  =  2 

0.0294 

0.0365 

0.0239 

0.0265 

0.0225 

Appendix  A 


IFTA  Design: 


\p{x,y)? 

x  =  -2 

X——1 

X  =  0 

X  =  1 

X  =  2 

y  =  -2 

0.0284 

0.0254 

0.0295 

0.0306 

0.0266 

y  =  -l 

0.0272 

0.0298 

0.0278 

0.0262 

0.0276 

y  =  0 

0.0254 

0.0296 

0.0289 

0.0296 

0.0254 

y  =  l 

0.0276 

0.0262 

0.0278 

0.0298 

0.0272 

y  =  2 

0.0266 

0.0306 

0.0295 

0.0254 

0.0284 

Annealed  Design: 


x  =  -2 

X  =  — 1 

X  =  0 

X  =  1 

x  =  2 

y  =  -2 

0.0316 

0.0315 

0.0331 

0.0285 

0.0285 

y  =  -i 

0.0302 

0.0318 

0.0297 

0.0297 

0.0313 

y  =  0 

0.0301 

0.0296 

0.0322 

0.0296 

0.0301 

y  =  l 

0.0313 

0.0297 

0.0297 

0.0318 

0.0302 

y  =  2 

0.0285 

0.0285 

0.0331 

0.0315 

0.0316 
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